UTHEP-317 
UTCCP-P-3 
August 1995 



QCD phase transition with two flavors of 
Wilson quarks using a RG improved action^ 



K. Kanaya 

Center for Computational Physics, Institute of Physics, 
University of Tsukuba, Ibaraki 305, Japan 

O 

o 

m 

^ Abstract 

^ The finite temperature QCD phase transition is studied on the lattice 

^ with degenerate two flavors of Wilson quarks. Motivated by reported strange 

O behaviors with the standard action on lattices with the temporal extention 

^ Ni = A and 6, a renormalization group improved gauge action is applied. On 

an A^i = 4 lattice, the strange behaviors observed with the standard action 
^ are removed with our improved action. The finite temperature transition 

is continuous in the chiral limit and it becomes quite smooth in all observ- 
I ables we studied when we increase the quark mass by increasing (5 along the 

^ crossover line. 
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1 Introduction and motivations 



Understanding the natiire of the finite temperature transition in QCD with 
two degenerate light quarks (Np = 2) is a very important step towards the 
clarification of the QCD transition in the real world. For Np = 2 QCD, 
Wilczek and Rajagopal have developed an argument that, if the transition in 
the chiral limit ("chiral transition") is of second order, then it will belong to 
the same univesality class as the three dimensional 0(4) Heisenberg model 
[1]. This gives us several useful scaling relations near the chiral transition 
that we can use to discriminate the nature of the chiral transition. If the 
chiral transition is of second order, we expect that it turns into an analytic 
crossover at finite quark mass. 

From numerical studies using staggered quarks [2] no signs of discontinu- 
ities are found for > 0.0125 on Nf = 4 and 6 lattices, where Nf is the 
temporal lattice size. Fiirthermore, recent scaling studies [3] show that be- 
haviors of physical quantities are consistent with the predicted 0(4) scaling 
near the chiral transition, suggesting that the chiral transition is of second 
order. However, because the action for Np = 2 staggered quarks is not local, 
further study with Wilson quarks is desired. 

The phase diagram for Np = 2 QCD with Wilson quarks is siimmarized 
in Fig. 1 collecting recent data of various groups (see references in [4]). As 
in the case of staggered quarks, the chiral transition at /3 = pet is shown 
to be continuous both on A^^ = 4 and 6 lattices [5, 6]. However, a detailed 
study by MILC collaboration [7] show that on A^^ = 4 lattices the transition 
becomes once very strong at A' ~ 0.19 when we decrease K along the finite 
temperature transition/crossover Kt, and it becomes weaker again when we 
further decrease K. Fig. 2 shows at /3 = 5.0 obtained by us and MILC 
collaboration. We find that changes its behavior sharply at K^. Other 
quantities, such as Polyakov loop and plaquette, also show a rapid change of 
behaviors at K^. Although no discontinuities are observed on A^ = 4 lattices, 
the transition does not look like a crossover here. On an A^ = 6 lattice [8] 
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MILC group even found first order signals at several intermediate values of 
K (K = 0.17, 0.18 and 0.19). Consulting the phase diagram Fig. 1 more 
closely, we, however, note that these points of strong are just the points 
where Kt lines get once very close to line due to the rather sharp bend 
of Kr at /3 5 caused by the cross-over phenomenon between the weak and 
strong coupling regions of QCD. Therefore, it seems plausible that the strong 
Kt is a lattice artifact caused by this unusual relation of Kt and K^ lines. 

Another lattice artifact we encounter with Wilson quarks using the stan- 
dard action is the strange behavior of quark mass rUg at /3 < 5.3 in the de- 
confined phase [9, 7]. for Wilson quark is defined through an axial Ward 
identity [10]. Numerical studies for /3 = 5.85 in quenched QCD [11] and 
P = 5.5 in Np = 2 QCD [12] show that (i) at zero temperature vanishes 
at K which is almost identical to that where = 0, and (ii) the value of 
mq does not depend on Nt and the phase. At /3<5.3, however, the second 
property is lost, in the deconfined phase now shows dependence on 
and deviate from the value in the confining phase (Fig. 2) suggesting sizable 
0(a) corrections to the axial Ward identity in the deconfined phase. 

A naive way out of these lattice artifacts is to increase so that we 
have the chiral transition in the weak coupling region, P^t > 5.5. However, 
our previous study on a 18^ x 24 lattice suggests that this requires A^^ > 18 
[6]. 

This motivated us to study these issues using improved actions. We apply 
a renormalization group (RG) improved action proposed by Iwasaki about 
ten years ago [13] which is reviewed in the next section. In section 3, results 
of numerical simulation performed on 8^ x 4 and 8^ lattices are summarized. 
We find that the unsatisfactory features of the standard action discussed 
above are removed at /3 down to the chiral transition point for Nt = 4. 



2 



2 Improved action 



In the scaling region, a lattice action describes the continuum properties at 
sufficiently large distances in lattice units. For these distances, the symme- 
tries expected in the continuum limt, such as the rotation symmetry, will 
recover. At short distances comparable with the lattice spacing a, however, 
lattice actions generally fail to reproduce continuum properties even in the 
scaling region. This makes MC simulations near the continuum limit very 
hard. On the other hand, lattice action has freedom to introduce less local 
terms without affecting the continuum limit. For example, lattice gauge ac- 
tion may contain rectangles, chairs, etc. besides the plaquettes adopted in 
the standard action: from loops up to length 6, 

S gauge — /3{Co X/ I 1 

where /3 = and Cq -|- 8ci -|- 16c2 + 8C3 = 1. Althoiigh the long distance 

behaviors are not affected by these additional terms, short distance properties 
are modified. Therefore, by adjusting these additional coupling parameters, 
the short distance lattice artifacts may be reduced. Such actions are called 
as improved actions. 

Several strategies have been proposed to find a set of coupling parameters 
which improves the theory efficiently. Here, we apply a strategy based on a 
renormalization group (RG) transformation [13]. 

Consider a block transformation of scale factor 2. Repeated applica- 
tion of a block transformation defines a chain of effective actions (trajectory) 
approaching an IR fixed point Pq in the infinite-dimensional coupling param- 
eter space (Fig. 3). The hypersurface where the correlation length diverges 
is called as "critical surface". We suppose a fixed point P^c on the critical 
surface. The trajectory connecting and Pg is called as "renormalized tra- 
jectory" (RT). Actions on RT are "perfect" (completely improved) actions 
[14] in the sense that continuum behavior is realized from the distance of 
la (up to distortions depending on the choice of the block transformation) 



3 



because RT is connected with the continuum theory P^c- When the initial 
action 5**^' is in the scaling region that is close to the critical surface, the 
trajectory first flows along the critical surface and then leaves the critical 
surface near Pqc towards Pq approaching RT gradually. The number of block 
transformations required to get sufficiently close to RT is nothing but log2 
of the minimum distance to obtain continuum behaviors with S^'-'K 

Hasenfratz and Niedermayer obtained a functional equation for Poo in 
several asymptotically free theories [14] and proposed to use it (approximated 
on a finite dimensional parameter space) as an approximate perfect action. 

Here, we are not aiming at a perfect action. Instead, we are trying to take 
advantage of the extended parameter space to accelerate the approach to RT. 
For example, consider an action Simp in fhe extended parameter space and 
suppose that the additional parameters in Si^p are adjusted so that the action 
gets sufficiently close to RT after, say, 2 block transformations. This means 
that, when we simulate the system using Simp, a separation of 4a is enough 
to get the continuum properties. In summary, our strategy is: (i) restrict 
ourselves to a small dimensional coupling parameter space consisting only 
with interactions which are easy to be implemented in numerical simulations, 
and (ii) find coupling parameters which minimizes the distance to RT after 
a few block transformations. 

Iwasaki applied this program to the SU(3) gauge theory [13]. Because 
the theory is an asymptotically free theory, correlation functions close to 
the critical surface can be computed by perturbation theory: for Wilson 
loops W{C) = 1 - g^F{C) + 0{g*), the first coeflacient F{C) can be easily 
computed. Choosing the block transformation sufficiently simple, blocked 
correlation functions can also be computed: W^^\C) = 1- g^F^^\C) + 0{g*) 
for J-th blocked lattice, even in the limit of RT (J = oo). Iwasaki studied 
the 4 parameter gauge action consisting loops up to length 6. To the lowest 
order perturbation theory C2 and C3 appear only in the combination C2 + C3. 
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Using a simple block averaging for the block transformation: 



where A^J\n) is the vector potential on /-th blocked lattice {n} and n' is a 
hypercube, and by defining the distance to RT after /-blockings by 



Ci < and C2 + C3 < 0, the valleys are given by 1.531ci + 02 + 03 = —0.5061 
and 1.761ci + C2 + C3 = —0.5116 for R^^^ and respectively, where the 

true minima locate at about Ci = —0.27 for both cases. Smallness of the 
resulting minimum values of R^^^ can be confirmed by the stability of F^^^'s 
under block transformations. 

The flatness of the bottom of these valleys suggests us to restrict the 
coupling parameter space further to C2 = C3 = 0. We then get a quite simple 
action consisting only with plaquettes and rectangles: 

Simp = P{Cq J2 I— I 

with Co = 1 - 8ci and ci = -0.331 (-0.293) to minimize i?^^) (i?^^)). In the 
foUowings we use Ci = —0.331. 

A perturbative calculation shows that the scale parameter of this action, 
AjM, is close to that in the MS scheme: Ajfg/ Ajm — 0.488 [15] in contrast 
with the large ratio for the case of the standard action. Therefore, no tadpole 
improvements [16] are necessary for our action. Correspondingly, ~ 1.8 - 
1.9 GeV is obtained at small /3 = 2.4 with the improved action, while /3 ~ 5.3 
is reqiiired with the standard action. This action was applied to mimerical 
studies of the string tension, the hadron spectrum, topological properties and 
the U{1) problem[17], and the resiilts indicate that 0(a) effects with this ac- 
tion are smaller than the standard action. Application to finite temperature 
gauge theory is in progress. 




with J2c over loops up to length 6, Iwasaki found that both R^^^ and i?^^^ 
form very narrow valley with rather flat bottom in the space (ci, C2 + C3): For 
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3 Results for A^^ = 2 



As a first step toward an improved full QCD, we combine the RG improved 
pure gauge action discussed in the previous section with the standard Wilson 
quark action. This should work well at least at small K where the effects 
of quarks can be absorbed by renormalizations of pure gauge interactions. 
Our results for Np = 2 QCD presented below suggests, however, that our 
action improves the theory very much also near Kc- Preliminary results are 
reported in [9, 4]. 

Simulations are done on 8^ x 4 and 8'^ lattices. Fig. 4 shows the result for 
the phase diagram. The chiral limit Kc is determined by = on the low 
temperature S"* lattice. Finite temperature transition line is determined 
on the 8^ X 4 lattice. Unlike for the standard action (Fig. 1), Kt goes simply 
away from K^, when we increase p. 

Shown in Fig. 5 is and rriq as a function of 1/K — l/K^. for various 
p. From the same a^^ ~ 1 GeV, results at /? 2.0 should be compared with 
those at /3 = 5 with the standard action shown in Fig. 2. The straight line 
envelop of in Fig. 5 corresponds to the chiral behavior oc rUq in the 
confining phase and deviation from this line signals the transition Kf to the 
deconfining phase. We note: (1) In contrast with the case of the standard 
action, shows no phase-dependence. (2) at decreases smoothly to 
zero when we decrease j3 toward P^t- This suggests that the chiral transition 
is continuous as in the case of the standard action. (3) At p > p^, as well 
as other physical observables are quite smooth at Kt. This suggests that the 
transition there is actually a crossover, in accord with a naive expectation 
that the second order chiral transition will become a crossover at finite m^. 

4 Conclusion 

Phase structure of finite temperature Np = 2 QCD with Wilson quarks is 
studied using a RG improved action. As in the case of the standard action. 
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the chiral transition is shown to be continiioiis. On the other hand, not Hke 
the case of the standard action where lattice artifacts make the transition very 
strong at intermediate tlie transition witli the improved action becomes 
quickly smooth when we decrease in accordance with a naive expectation 
about the fate of a second order chiral transition at finite rUq. Studies for 
Np = 3 and 2 + 1 are in progress. 

Simulations are performed with HITAC S820/80 at KEK and with QCD- 
PAX and VPP500/30 at the University of Tsukuba. I would like to thank 
members of KEK and the members of QCDPAX collaboration for their hospi- 
tality and strong support. This work is in part supported by the Grant-in- Aid 
of Ministry of Education, Science and Culture (No.07NP0401). 
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Figure 1: Phase diagram for A'^^ = 2 QCD with Wilson quarks using the 
standard action. K^. is the chiral hmit on T = lattices and Kt(Ntys are 
for finite temperature transition/crossover. /5rf(A^t)'s are for crossing points 
of Kt and Kc lines, i.e. the chiral transition. Lines are for guiding eyes. 
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Figure 2: and 2mq with the standard action at /3 = 5.0. 
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Figure 3: Renormalization group flows in the infinite-dimensional coupling 
parameter space of QCD. 
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0.24 




Figure 4: Phase diagram for Np = 2 QCD with Wilson quarks coupled to an 
RG improved gauge action. is the chiral limit determined by = on 
an 8'^ lattice and Kf is the finite temperature transition/crossover obtained 
on 8^ X 4. Lines are for guiding eyes. 
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gure 5: and 2mg with a RG improved action at various 
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